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ABSTRACT 

Recent discovery of the year-scale variability in the synchrotron X-ray emission of the supernova remnant 
(SNR) RX J17 13.7-3946 has initiated our study of multi-epoch X-ray images and spectra of the young SNR 
Cassiopeia A based on the Chandra archive data taken in 2000, 2002, and 2004. We have found year-scale time 
variations in the X-ray intensity for a number of X-ray filaments or knots associated with the reverse- shocked 
regions. The X-ray spectra of the variable filaments are characterized by a featureless continuum, and described 
by a power law with a photon index within 1.9-2.3. The upper limits on the iron K-line equivalent width are 
110 eV, which favors a synchrotron origin of the X-ray emission. The characteristic variability timescale of 4 
yr can be explained by the effects of fast synchrotron cooling and diffusive shock acceleration with a plausible 
magnetic field of 1 mG. The X-ray variability provides a new effective way of studying particle acceleration at 
supernova shocks. 

Subject headings: acceleration of particles — ISM: individual(Cassiopeia A) — radiation mechanisms: non- 
thermal 



1. INTRODUCTION 

The young (~ 330 yr old) supernova remnant (SNR) Cas- 
siopeia A (Cas A) is a unique astrophysical laboratory for 
studying high-energy phenomena in SNRs because of its 
brightness across the entire e lectromagnetic spectrum from 
radio (e.g.. iBaars et al.1 Il977 | ) to verv-hi^h-e nergy gamma- 
rays (I Aharonian et al.ll20Qlb lAlbert et alJl2007b , While the X- 
ray emission of Cas A below 3 keV is dominated by line emis- 
sions from the thin thermal plasma of the shocked ejecta, as 
demonstrated by observations with ASCA (Ho lt et al"1ll994|) . a 
prominent nonther mal tail is known t o be present and d omi- 
nant above 10 keV (lAUen et al.lll997UFavata et al.lll997l) . 

Nonthermal X-rays appear to be emitted mainly from frag- 
mented filaments or knots in reverse- shocked regions. The 
bulk of the hard X-ray emission in 10-12 keV from Cas A 
comes from such components; the total e mission in the outer 
peripheral filaments f ound by Chandra (Hug hes et al.l 2000: 
IVink & La ming 2003|) i s only a small fracti on of the total 4- 
6 keV continuum flux (Bleeker et al. 2001). The interpreta- 
tion of the nonthermal com ponent in reve rse- shocked ejecta 
remains an open issue (e.g.. lLaming|l200 l[) . 

In this Letter, we present the analysis of multi-epoch Chan- 
dra imaging and spectrometric data of Cas A. This work 
was motivated by our recent discovery of time variability in 
the synchrotron X-ray emission of SNR RX J1713.7-3946 
(lUchiyama et al.ll2007|) . where the year- scale time variability 
requires, most likely, the presence of a largely amplified mag- 
netic field of 1 mG in compact filaments. The strength of mag- 
netic field in t he compact radio k nots in Cas A is estimated 
to be 1-3 mG (lAtoyan et al.ll2000|) . just sufficient to produce 
year- scale variability through synchrotron cooling and shock 
acceleration of multi-TeV electrons. Indeed, here we report 
significant time variations in nonthermal X-ray emission from 
many compact regions. The distance to Cas A is assumed to 
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TABLE 1 

LOG OF Chandra OBSERVATIONS 



ID 


OBSID 


Date 


R.A. 


Decl. 


Roll 


Exposure 








(deg) 


(deg) 


(deg) 


(ks) 


PI 


114 


2000 01 30 


350.9176 


58.7934 


323.4 


49.9 


P2 


1952 


2002 02 06 


350.9172 


58.7938 


323.4 


49.6 


P3 


5196 


2004 02 08 


350.9188 


58.7949 


325.5 


49.5 



NOTE. — Sky coordinates are right ascension (R.A.) and decli- 
nation (Decl.) in J2000.0, representing the pointing direction. The 
exposure time is the net integration time after standard filtering. 

be 3.4 kpc dReed et al.ll 19951) so the angular scale of V corre- 
sponds to 1 pc. 

2. DATA 

We analyzed the X-ray data obtained with the ACIS-S3 
chip on 2000 January 30 (Hwang et al. 2000; Gotthel fetaD 
2001), 2002 February 6 foeLanev & Rudnickl l2003h . and 
2004 February 8 dHwang et al .112004 . each with ~ 50 ksec 
exposure. Table [T] gives the log of the Chandra observations. 
The observations were made with ACIS-S3 in the graded 
mode, with which telemetry size can be reduced. The three 
observations have almost identical configuration such as the 
aim point and the roll angle, allowing us to directly com- 
pare the three images without taking account of the varying 
point- spread functions (PSFs). The aim point is located ~ 1 .9' 
southeast of the central point source, the compact remnant left 
behind after the supernova explosion. 

The X-ray analysis was done using the CIAO software 
package (version 3.4). We reprocessed the event data, ap- 
plying standard data reduction with the latest calibration files. 
We find that the sky coordinates of the central point source 
in the three epochs are consistent with each other with an ac- 
curacy of < 0.2". The point source itself is estimated to be 
moving at an angular speed of ~ 0."02 yr -1 dThorstensen et al.l 
12001b . thus 0."08 in the 4-yr span. 

Figure Q] shows the Chandra image (with 0."2 x 0."2 pix- 
els) of Cas A in the 4-6 keV band averaged over the three 
epochs. This energy band corresponds to the continuum 
emission between Ca and Fe K lines, comprised of thermal 
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Fig. 1 . — Chandra ACIS-S image of Cas A in the 4-6 keV band av- 
eraged over the three epochs. North is up and east is to the left. 
The image is displayed with a square root scaling in a range of 
0-3 counts pix -1 . The outer circle with 153" radius denotes the for- 
ward shock positions, and the inner cir cle with 95" radius i ndicates 
the mean location of the reverse shock dGotthelf et al .11200 II) . 

bremsstrahlung and nonthermal component(s). Assuming that 
the hard X-ray above 10 keV is dominated by nonthermal ra- 
diation processes, we estimate about 50% of 4-6 keV photons 
should be of nonthermal origin. Filamentary structures hav- 
ing high surface brightness in the continuum band are found 
only for radii less than 2' (referred as inner filaments), imply- 
ing their as sociation with the rev erse shock as is shown to be 
the case bv lHelder & Vinkl (120081) . 

We are interested here in the nonthermal X-ray compo- 
nents. We selected the two boxes shown in Fig. Q] based on 
the following criteria. First, we identified regions sufficiently 
bright in the 4-6 keV band, exceeding 3 counts pix -1 at least 
in one epoch, in order to obtain statistics enough for a variabil- 
ity search. Note that this condition is met only by the inner 
filaments, associated with the reverse shock or its secondary 
shocks. The second requirement is that the level of the con- 
tinuum emission is high relative to the silicon line strength. 
Specifically, we calculated the ratio of 4-6 keV counts to 1.7- 
2.2 keV counts at each pixel, and selected those regions hav- 
ing the ratio exceeding 0.5. In this way, we found that the 
inner filaments satisfying these criteria are distributed inside 
the two boxes only. Since we do not aim at constructing a sta- 
tistical sample in this work, we adopted these simple criteria, 
accounting for neither position-dependent PSFs nor effective 
area. 

3. ANALYSIS 

Figure [2] displays a sequence of three-epoch 4-6 keV im- 
ages of the central box (left) and the western box (right), 
each with dimensions of 0.5' x V . Morphologically distinc- 
tive components that satisfy the criteria mentioned above are 
labeled as A, B, C, E, G, and H. Also we tagged an additional 
two features, D and F. 

The 4-6 keV images over a 4-yr time interval reveal a dozen 
time- variable structures. Region C has decreased by a fac- 
tor of 2 from 2000 to 2004. Region H underwent a flux 
brightening by 50% over the time interval of 4 yr. (By mea- 
suring the silicon line flux from the compact ejecta features, 



TABLE 2 
Results of Spectral Fitting 



ID 




r 


^2-10 


X 2 (d.o.f.) 




10 22 cm" 2 




10 _12 erg cm -2 s _1 


AB 


2.00±0.10 


2.31 ±0.08 


1.59 ±0.04 


0.85 (274) 


C 


1.80±0.13 


2.37±0.11 


0.77 ±0.03 


1.08 (221) 


DEFG 


1.69 ±0.06 


2.24 ±0.05 


5.21 ±0.05 


0.98 (365) 


H a 


1.80 (fix) 


1.97±0.10 


3.06 ±0.07 


0.84 (206) 



NOTE. — Fitting X-ray spectrum by a power law with photoelec- 
tric absorption. Absorbing column density Mi, photon index T, and 
unabsorbed 2-10 keV flux F2-10 are shown with 90% errors. 
a X-ray energy below 3 keV is excluded from the fit and the absorb- 
ing column density is fixed at Mi = 1.8 x 10 22 cm -2 . 



we estimated a systematic error in the flux measurements to 
be < 9%.) The variable filaments are spatially extended be- 
yond Chanda's PSFs, typically ~ 3" -5" in length and ~ I" 
in width. Also, some filaments that we did not label, such as 
those in the south of region C, exhibit time variability. Year- 
scale time variability appears to be the prevailing nature in the 
bright filaments or knots in the 4-6 keV band. 

Before investigating the epoch-dependence of the spectral 
parameters, we first constructed X-ray spectra averaged over 
the three epochs. Regions A and B were combined into a 
single zone, and regions D, E, F, and G were merged. The 
background spectrum taken from a nearby, low brightness re- 
gion was subtracted from each spectrum. Since the choice 
of a background region largely affects the low-energy part of 
the X-ray spectrum of region H, photons below 3 keV were 
discarded for this region. The four X-ray spectra (A+B, C, 
D+E+F+G, and H) in the 1-8 keV band are shown in Figured 
They all appear to be line-less. 

To characterize the shape of the line-less X-ray spectrum, 
we adopted a simple power-law model attenuated by inter- 
stellar absorption, which yielded statistically good fits. We 
obtained the best-fit parameters in ranges of Y ~ 2.0-2 .4 and 
N H = 1.7-2.0 x 10 22 cm" 2 as summarized in Table [2 The 
absorbing column densit y is broadly consistent with those 
for thermal components ( Willingale et al.l 120021) . We also 
performed spectral fitting using a thermal bremsstrahlung 
model. A statistically acceptable fit was obtained in each 
case. The best-fit electron temperature is in a range of 4.0-6.6 
keV, much higher than the electron temperatures, 1-2 keV, 
derived for the X-ray k nots showing thermal spectra (e.g., 
lLaming & Hw ang 2003). Moreover, by combining all the 
spectral data, the equivalent width of the iron K lines is con- 
strained to be < 110 eV. The absence of the emission lines dis- 
favors any bremsstrahlung interpretation of the featureless X- 
ray s pectra, includi ng bremsstrahlung by suprathermal elec- 
trons (lLamingll200l|) . 

In order to track the time variations in terms of both inten- 
sity and spectral shape, the X-ray spectrum extracted from 
each labeled region at each epoch was fit by an absorbed 
power-law model. The absorbing column density was fixed at 
a relevant value in Table [2 Figured shows the time sequence 
of the unabsorbed X-ray flux normalized by that in the year 
2000. Four out of the eight regions exhibit large time varia- 
tions (> 10% yr -1 ). The time sequence of photon index Y for 
each region is also presented in Figure 01 We found notable 
variations of Y only in two cases, regions F and H. Region H 
gives the hardest index of Y = 1.85 ± 0.10 in 2000, suggest- 
ing a direct link between flux increase and spectral hardening. 
The change in Y found in region F should be treated with care, 



Synchrotron X-ray Variability in Cas A 
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Fig. 2. — A sequence of three-epoch 4-6 keV images of the two 0.5' x 1' boxes in Fig.[T] The images are shown in a linear scale in a range 
of 0-3 counts pix -1 for the left panel, and 0-4 counts pix -1 for the right panel, respectively. Pixels have dimensions of 0."2 x 0"2. Gaussian 
smoothing with a kernel of 0"8 is applied. The central box (i.e., the left panel) is close to the aim point and therefore the PSF is sharp as 
is evident from the point source. The western box (the right panel) is away from the aim point, so that some of the spatial extent should be 
attributed to the broadening by the PSF. 



> 

M 




Fig 



channel energy (keV) 
3. — X-ray spectra of combined regions: AB (black), C (red), 



DEFG (green), and H (blue). The three epochs are averaged. The 
best-fit power-law function with interstellar attenuation is superim- 
posed in each case. Nearby diffuse X-rays are subtracted from each 
source spectrum. 

because region F faded after 2000 and became almost invis- 
ible in 2002 and 2004. Photon index measured in 2002 and 
2004 would not be physically connected with that in 2000. 

4. DISCUSSION 

We have found that a half (4/8) of the filamentary regions 
that were selected simply by their brightness and hardness, 
showed a year- scale variability in X-ray flux. The X-ray spec- 
tra of the selected regions in the 1-8 keV range are featureless 
and well fit by an absorbed power law with a photon index 
T ~ 1.9-2.4. Given the absence of emission lines, we pre- 
fer to attribute the nonthermal X-ray emission in regions A-H 
to synchrotron radiation by very high-energy electrons. Syn- 
chrotron X-ray filaments or knots would originate in slightly 
overdense ejecta, say by a factor of ~ 2, into which a reverse 
shock (or transmitted shocks) with a typical shock speed o f 
^ 1000 km s -1 has been driven ( see lLaming & Hwa ng 2003). 

Recently, Patnaude & Fesen d2007h have reported X-ray 
variability in compact thermal knots, resulting from the pas- 
sage of reverse shock over the ejecta clumps. In the entire face 
of the remnant, they identified six time- varying structures; 
all but one region (their R4, corresponding to our region C) 
were modeled by a "thermal" spectrum, specifically by a non- 
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Fig. 4. — Time sequence of the unabsorbed 2-10 keV flux relative 
to that in the year 2000 for eight filamentary regions (filled circles), 
as well as time sequence of photon index V (open circles) are shown 
with la errors. 

equilibrium-ionization plasma with an electron temperature 
of ~ 1 keV. On the other hand, we found a comparable num- 
ber of variable "nonthermal" features just inside the two small 
boxes, by making use of 4-6 keV maps which are sensitive to 
nonthermal variability. Morevoer, if we a pply our criteria of 
significant variability, > 10% yr -1 , to the Patnaud e~& Fesenl 
(2007) sample, only two features remain. We conclude that 
variable nonthermal features are much more prevailing than 
the thermal ones. 

Cas A is a very bright nonthermal radio source which indi- 
cates the presence in the remnant of relativistic electrons with 
a Lorentz factor 7 e > 100. The same electrons produce also 
high energy gamma-rays through relativistic bremsstrahlung. 
Since the gas density in Cas A is known quite well, the 
gamma-ray flux is robustly controlled by the strength of the 
magnetic field. Therefore, high energy gamma-ray observa- 
tions, combined with radio dat a, can provide reliable mea- 
surements of the magnetic field (ICowsik & Sarkar Ill980|) . A 
simple homogeneous model based on the radio observations 
and the upper limits of the > 100 MeV flux reported by the 
EGRET collaboration (Esp osito et al.l Il996|) yields a lower 
limit on the magnetic field of 0.4 mG ( Atovan et al. 2000). In 
a more sophisticated two-zone model of lAtoyan et al.l (|2000), 
the magnetic field in the compact radio knots, where particle 
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acceleration is supposed to be taking place, is constrained to 
be B = 1-3 mG. These values would reasonably be applicable 
to the X-ray filaments, which have a similar size with the ra- 
dio features, and are also presumed to be the sites of efficient 
particle acceleration by reverse (or secondary) shocks. 

It is interesting to ask whether the X-ray flux changes are 
attributable solely to the changes in magnetic field strength. 
If so, it is expected that a substantial flux change in syn- 
chrotron X-rays accompanies a marked change in radio inten- 
sity. While the radio flux is proportional to B 2 , the X-ray flux 
has weaker dependence on B, because X-rays are emitted in a 
quasi- saturated regime. Therefore one would expect stronger 
va riations in the radio t h an in the X-ray band. According 
to lAnderson & Rudnickl (|1995|) . however, among about 70 
bright radio compact knots exceeding 20 mJy beam -1 , none 
showed a brightness increase larger than +10% yr -1 in the 
archive Very Large Array data at 5 GHz. A typical rate of 
brightening is about +3% yr -1 , corresponding to an ^-folding 
time of ~ 30 yr, which is much longer than the X-ray vari- 
ability timescale of 4 yr. Therefore it is more likely that the 
observed time variations are caused mainly by a rapid increase 
or decrease in the density of X-ray-emitting electrons. (The 
increase in the radio-electron density above the "sea" level is 
expected to be less significant, provided that fresh electrons 
have a harder energy spectrum than old electrons.) Indeed, 
as we argue below, the known constraint on B provides the 
synchrotron cooling time being compatible with the observed 
time variability. 

We have reported a similar X-ray variability in SNR 
RX J17 13.7-3946, and attributed it to fast s ynchrotron cool- 
ing a nd shock acceleration of electrons (lUchiyama et al.l 
2007). Let us examine the obser ved varia bility timescale, us- 
ing the equations presented in lUchiyama et al.l (|2007|) . Syn- 
chrotron cooling time is given by t syn ~ 1 .5 5~q 2 'e^v Y r ' 
where B m Q is the magnetic field strength in units of mG and 
£kev is the mean synchrotron photon energy in units of keV. 
For B ~ 0.5 mG, the synchrotron cooling timescale becomes 
comparable to the variability timescale, t 



syn 



4 yr. 

Therefore, even a magnetic field of 0.4 mG is able to explain 
the observed variability. For B ~2 mG, if the injection of rel- 
ativistic electrons responsible for the synchrotron X-rays has 
suddenly ceased in a filament, the X-ray flux can be reduced 
on a timescale of t syn ~ 0.5 yr. If the electron injection rate 
(i.e., the number of particles picked up and injected into the 
acceleration process in unit time) decreases noticeably in 4 yr, 



the number of X-ray-emitting electrons and consequently X- 
ray intensity can drop as observed. We note that plasma char- 
acteristics may change on a timescale of a few years given the 
arcsecond-scale size of the thermal ejecta filaments. 

Rising flux is a signature of ongoing production (accel- 
eration) of synchrotron-emitting electrons with energies of 
~ 10 TeV. It implies that electron energy can increase to 
~ 10 TeV in 4 yr or less at the shock. Assuming that elec- 
trons are accelerated via diffusive shock acceleration, the ac- 
celeration time of the X-ray-emitting electrons is given by 

face ~ 9 ^ 5 m 3 G 2£ k^v^foooy r ' where rj > 1 is a gyrofactor, and 
Viooo is the velocity of the reverse shock relative to the un- 
shocked ejecta ahead of it in units of 1000 km s" 1 . There- 
fore, reverse shocks with a typical speed of ~ 1000 km s -1 
(lLaming & Hw ang 2003) and a magnetic field of 2 mG can 
accelerate X-ray-emitting electrons within a timescale of 
^var ~ 4 yr, provided that the acceleration proceeds close to 
the Bohm diffusion limit (77 ~ 1). 

The presence of time variability together with the absence 
of the Fe-K lines strongly prefer a synchrotron origin of the 
hard X-ray emission in the selected filaments, providing new 
evidence of acceleration to multi-TeV energies at the reverse 
(and transmitted) shocks in ejecta. Though the regions we 
studied here contain only a portion of the reverse- shocked 
gas, the bulk of the hard X-rays above 10 keV is also likely 
of a synchrotron origin. If so, the steepness of the c ontin- 
uum spectrum above 10 keV, T ~ 3 dFavata et al.lll997|) . indi- 
cates that the synchrotron spectrum rolls off at around 1 keV, 
being con sistent with the theory of diffusive shock acceler- 
ation (see lUchiyama et al. 2007). The maximum energy of 
the electrons is thus constrained to be ~ 10 TeV. Therefore 
the TeV gamma-rays measured by HEGRA (lAharonian et al.l 
12001b and MAGIC d Albert et al. 1120071) groups are difficult 
to be explained by electron bremsstrahlung in the reverse- 
shocked ejecta, and are likely to be of a hadronic origin. 

The interpretation of the X-ray variability has to be crit- 
ically assessed through future theoretical and observational 
works. Regardless of what is behind the nonthermal X-ray 
variability, it offers us a new diagnostic tool to understand 
particle acceleration in supernova shocks. 



We wish to thank Herman Lee for his useful comments on 
the original manuscript. 
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